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Abstract 

We study a mechanism for superconductivity in quasi-one-dimensional ma- 
terials with Ising anisotropy. In an isolated chain Ising anisotropy opens a 
spin gap; if inter-chain coupling is sufficiently weak, single particle hopping is 
suppressed and the physics of coupled chains is controlled by a competition 
between pair hopping and exchange interaction. Spin density wave and triplet 
superconductivity phases are found separated by a first order phase transi- 
tion. For particular parameter values a second order transition described by 
SO (4) symmetry is found. 



I. INTRODUCTION 

This paper study a mechanism of superconductivity for quasi-one-dimensional (Q1D) 
metals. Well known examples of Q1D materials exhibiting superconductivity are organic 
compounds (TMTSF^X and (TMTTF^X. At this moment there is no consensus on the 
origin of the superconductivity this metals. Two possibilities are discussed in the litera- 
ture. One is BCS theory with the attraction mediated by either phonons or magnons |T],@. 
The second approach developed by early workers in the field tries to accommodate the 'g- 
ology' picture of a one-dimensional conductor to the materials of interest ||. The crystal 
structure of the materials allows one to treat them as an array of the weakly coupled one- 
dimensional conductors formed by chains of flat molecules. It is assumed further that at 
hight temperature every conductor can be approximately described as a Luttinger liquid 
(LL) characterized by dimensionless parameters /C CjS . Depending on those parameters, at 
least one of four susceptibilities, spin-density, charge- density, singlet Cooper pair or triplet 
Cooper pair, diverges at low temperature. In a three-dimensional crystal a phase with a 
broken symmetry corresponding to the most divergent susceptibility will be formed at low 
temperature. Thus, the efforts were invested into the study of the LL parameters as a func- 
tion of the temperature and the pressure with the hope of finding conditions under which 
the Cooper pair susceptibility becomes divergent. Due to Coulomb repulsion, however, it 
was proven to be difficult to move the system into the region of divergent Cooper pair 
susceptibility HI. 
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The inability of the 'g-ology' picture to offer an explanation for the superconducting 
phase, by contrast to its success in describing the spin-density wave (SDW) phase |||| 
prompted many researchers to think in terms of a hybrid approach |]J. It is known that in 
the majority of Q1D materials SDW and Pierls-SDW phases are stable at ambient pressure. 
In order to stabilize superconductivity external pressure has to be applied. The pressure 
acts to increase the transverse coupling between chains driving the system closer to an 
anisotropic 3D Fermi liquid. It was suggested that this 3D Fermi liquid is a crucial require- 
ment for the superconductivity while the SDW origin is essentially one-dimensional. With 
this agreed upon one can use BCS theory and its strong-coupling modification to derive 
different experimentally relevant quantities. 

In this paper we will not follow this path. Instead, it will be demonstrated that both 
SDW and superconductivity can be derived within the framework of weakly coupled one- 
dimensional conductors. Our proposal is very similar in spirit, although different in detail 
to that of Carlson et. al. who studied consequences of stripe formation in models of the 
cuprate superconductors. 

In its current state the mechanism discussed in this paper cannot be applied directly to 
the organic superconductors because it requires substantial Ising anisotropy whereas in the 
organic materials deviations from the rotational symmetry are very small. Yet, it may be 
a useful step toward constructing viable mechanism for this materials. We discuss possible 
modifications at the end of the paper. It is also can be considered as a guide in the search for 
new superconducting materials. In addition to that it is interesting from purely theoretical 
point of view. 

There are two parts to the proposed mechanism. At high energies inter-chain couplings 
can be neglected and the physics is purely one-dimensional. Unlike the 'g-ology' picture, 
however, it is accepted that chains are not simple Luttinger liquids in high energy regime. 
Instead, we assume that an interaction which is anisotropic in spin space opens a gap A s in 
the spin sector. 

Opening of the gap dramatically changes the low energy behavior of the system. For 
reasonable values of /C c (1/2 < JC C < 1) there are now two relevant operators ||. One 
is longitudinal SDW (SDW Z ), the other is triplet superconductivity (TSC) dual to SDW. 
Provided that inter-chain transversal coupling is smaller then A s an effective low energy 
hamiltonian can be easily derived. It describes the competition of this two order parameters. 
The intra-chain electron-electron repulsion makes JC C smaller then unity. In this situation 
SDW is more relevant then TSC. Yet, under pressure, the influence of the next-to-nearest- 
neighbor coupling increases. This increases the effective Josephson coupling between the 
chains and decreases the effective exchange coupling. Therefore, it is possible to observe 
superconductivity under pressure. 

An essential feature of the approach is its many-body character. Its main ingredient is 
severe spin-charge separation. Unlike the usual BCS mechanism, this route to superconduc- 
tivity does not require attractive interaction. 

Another interesting property of the mechanism is that for some parameter values the 
symmetry group of the low energy effective hamiltonian is enhanced. Generically, the sym- 
metry is U(l)(gU(l), one copy of U(l) for the gauge transformation and another one for 
translation along the chains. If parameters of the hamiltonian are fine-tuned the symmetry 
group is enlarged to SO (4). An extra symmetry is particle-hole rotation connecting SDW 
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and TSC ground states. 

The rest of the paper is organized as follows. In Section II we introduce high energy 
single-chain hamiltonian and discuss its spectrum and correlation functions. We introduce 
our low-energy hamiltonian in Section III. The phase diagram is calculated in Section IV. 
The effect of long range Coulomb interaction is considered in Section V. The SO(4) symmetry 
is derived in Section VI. Finally, Section VII contains discussion. 



II. HIGH ENERGY DESCRIPTION OF THE SYSTEM 

It is usually agreed upon that high energy physics of Q1D materials is purely one- 
dimensional. The hamiltonian of a single chain can be written as follows: 

H = J L dxH, (1) 

n = n c + n s (2) 



n c = ^ (ic- 1 (v$ c ) 2 + jc a (ve c ) 2 ) (3) 

H s = I (JC; 1 (V$ s ) 2 + K B (V0 S ) 2 ) + ^ Vtmmim cos v/M s , (4) 
where we used the Abelian bosonization prescription B: 



to express the electron hamiltonian in terms of bosonic fields. In the above formulae f] pa 
are Klein factors, 9 S C are spin and charge boson fields, $ SjC are dual spin and charge fields. 
The index p = ± characterizes the electron direction of motion, '+' corresponds to 'L(eft)'- 
moving electrons and '-' corresponds to 'R(ight) '-moving electrons. At high energy the spin 
and charge bosons are decoupled. The parameters of the hamiltonian (|J) are phenomeno- 
logical quantities, functions of the microscopic electron interactions and dispersion. 

Equation (|3|) implies that the spectrum of charge excitations is gapless. The presence or 
the absence of the gap in the spin sector is a subtler matter, and depends on the relation 
between g bs and /C s . A renormalization group (RG) analysis of (|) leads to P,p!D|: 



= -2g b , (K, - 1) . (7) 
The infrared properties in the spin sector are defined by the value of the RG flow invariant: 

I = A~ ~ !) 2 ■ (8) 

For positive I the cosine term is relevant and the spin bosons spectrum has a gap. For I < 
the cosine term scales down to zero and the spectrum is gapless. Physically, the value of 
/ corresponds to the spin rotational symmetry of the interaction. If the interaction is easy 
axis (plane) then / > (J < 0). 

Let us now assume that the interaction has an easy axis. In this case there is a gap 
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A s ocAexp (-V 1/2 ) 



(9) 



in the spin sector. Here A is the ultraviolet cut-off energy. 

Opening of the spin gap does not imply long-range ordering. All expectation values 
disallowed by symmetries of the hamiltonian (p]) remain zero. For example, the operator 
density corresponding to SDW order takes the form: 

5,(2**) = i (>i T ^ RT - ^VrJ oc e 1 ^ sin v& s . (10) 

The non- vanishing spin gap fixes the field a/27t$ s near ±7r/2 so the expectation value 
(sin \/2ir& s ) is non-zero. The SDW order parameter, however, is zero because charge bo- 
son field $ c fluctuates: (S z ) oc (exp(i\/27r<I> c )) = 0. The correlation function (S Z S Z ) oc 
(exp(iv / 27r$ c (x, t)) exp(— i^/2n& c (x' , r'))) decays algebraically with the exponent /C c for 
large separations. 

The state with the spin gap was proposed and described first in Ref. ||. 

III. INTERCHAIN COUPLING FOR u < A s 

In this section we determine the effect of inter-chain coupling on the assumption uj < A s . 
The hamiltonian density which incorporates such interactions is: 

w = E^i+E w -L« ( n ) 

i ij 

where i and j labels chains and the single chain hamiltonian density (j|D is TCi. The single 
particle inter-chain hopping is represented by the hamiltonian density 7i±ij. It will be 
demonstrated below that the other physically important interaction the inter-chain density- 
density repulsion is irrelevant. 

An addition of a single electron or hole to a chain with a spin gap creates a state with 
a spin boson soliton. The energy of such soliton is at least A s . Yet, it is possible to add a 
pair of particles or a particle and a hole in such a way that the soliton is not created and 
the spin bosons are left undisturbed. This implies that while the single particle hopping 
in our situation is greatly diminished, the particle-particle (Cooper pair) hopping and the 
particle- hole (exchange) hopping survive the presence of the spin gap |7|,pTl|. Technically, 
when transverse hopping is smaller then A s the fast modes (those, whose energy exceeds 
the spin gap) can be easily 'integrated out'. The slow degrees of freedom are charge bosons 
with the energy less then the spin gap. The effective dynamics of these modes are given by 
the following hamiltonian density: 

K ff = |E (/C 1 (v$ cj ) 2 + k c (ve CJ ) 2 ) (is) 
z j 

+ 2a c' E ( J ?t (vUvk-iVk-iVi-i + '/[. . '/h , '//{ '//..) cos V2^ ($ ci - $ci) 

u 

- J!j (vUvkivkivli + vhvkAvU) cos V2^(® ci - Qcjj) ■ 
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Here J is the effective exchange coupling and J sc is the effective Josephson coupling. The 
hamiltonian itself is derived in the Appendix. The ultraviolet cut-off for this effective theory 
is equal to a^ 1 = A s /v c . 

It is necessary to explain the physical meaning of different terms of C[3~3|) . The origin 
of the first term is obvious. This term is responsible for the intra-chain dynamics of the 
charge bosons. The second term acts to order $ c . As one can infer from ([U]) the field $ c 
should be viewed as the phase of the SDW. When exp(i\/27r<I> c ) acquires finite expectation 
value the ground state becomes SDW^. One has to remember that (sm\/2iiQ s ) is non-zero. 
Unlike S z , the superconducting order parameter in SDW phase is zero. It is proportional to 
cxp Since G c is dual to $ c the ordering of the latter implies strong fluctuations 

of G c and, therefore, (exp (i\/27r©c)) = 0. 

The last term of ([13]) describes the Josephson coupling between pairs of the chains. If 
the field 6 C is ordered the ground state becomes the longitudinal triplet superconductor. 
The operator density 

A^o = VVil - ^kAi <* eiV ^ &C sin ( 14 ) 

acquires anomalous non-zero expectation value. We see that the field G c plays a role of the 
superconducting order parameter phase. The SDW order parameter has zero expectation 
value in the superconducting phase. 

In a purely one-dimensional system a ground state with a broken symmetry cannot exist. 
However, due to inter-chain coupling in ([13|) the problem at hand becomes three-dimensional 
where the broken symmetry phase can be stabilized at low temperature. In our situation 
the longitudinal triplet superconductivity and SDW 2 are the only types of broken symmetry 
which are allowed to exist. Other possibilities, such as charge- density wave, transversal 
SDW (SDW^,,), longitudinal TSC or singlet superconductivity are incompatible with the 
spin gap and (§). Which of the two allowed ground states, TSC or SDW, has lower energy 
at T = depends on the parameters of the effective hamiltonian. 

There is a relative minus sign between the second and the third term of (fl3|) . This 
has important consequence. If the ground state is a superconductor the order parameter 
is uniform over the chain array. In case of SDW, however, the order parameter has the 
opposite sign on the neighboring chains. Therefore, next to nearest neighbor coupling acts 
to stabilize the superconductivity and to destabilize SDW. This circumstance will be used 
when we will map out the phase diagram of ( |T3"D in the next section. 

Before we continue with the phase diagram we want to explain why the inter-chain 
density-density repulsion can be neglected. We have mentioned above that the charge- 
density expectation value is zero due to the spin gap. Therefore, in the first order of the 
inter-chain density- density coupling V± the contribution of this term is exactly zero. In the 
next order the contribution is non-zero. It is proportional to exp ^±i\/87r ($ C j — $ c j))- This 
operator is irrelevant and can be neglected. 

IV. PHASE DIAGRAM 

In this section we will determine the phase diagram of the system as a function of J's and 
temperature. Before we start constructing the diagram let us discuss possible constraints on 
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the values of the model parameters. There are several parameters our hamiltonian ([13]) is 
characterized by: a c , /C c , J^-. We assume that the exchange constants J sdw are equal to the 
Josephson coupling constants J sc . Why this is a reasonable approximation is discussed in 
the Appendix. We will put all equal to zero unless i and j are nearest neighbors or next- 
to-nearest neighbors. It is accepted that the next-to-nearest neighbor coupling constant J 2 
is smaller then the nearest neighbor constant J\\ J\> J 2 . When the interaction is repulsive 
/C c is bigger then one half. For the Hubbard model it has been found that /C c is smaller then 
two. We will confine /C c to this interval: 

1/2 < /C c < 1. (15) 

Outside of this interval the proposed mechanism cannot work: for /C c < 1/2 the Josephson 
coupling is irrelevant. 

To map out the phase diagram we will use mean field theory One constructs the mean- 
field hamiltonian for the coupled chain problem by re-writing the inter-chain interaction 
term: 

4 J2 Jij^ 1 cos v 7 ^ (©ci - O ci ) = Hmf + AH, (16) 
2 



'Hmf = ^2 "A? ( cos v^rOcj (cos \f2~nQ C j) + cos V27rQ C j (cos v^rO^ )J , (17) 



2 

AH = J ij ( 2 cos ^2tt (0 ci - Oci) ( 18 ) 



— cos \/27rO C j (cos \f2~K® c ^j — cos \f2~nQ C j (^cos v^tOc^) 
and neglecting ATi. The mean-field single-chain hamiltonian density is: 

HMF = ^2 H-ci ~ J-dc 1 COS V^TT&ci ~ J+O'c 1 COS V^TrQci, (19) 

% 

j_ = J_(cosV27r® c ), (20) 

J + = J + (cosV27r© c ). (21) 

The expression for 7i c is given by fl3|). Both cosines are relevant interactions. Therefore, 
they produce gap A c in the excitation spectrum of the charge boson. But the nature of the 
ground state, superconducting or magnetic, depends on the values of fC c and the effective 
coupling constants ('+' for TSC and '-' for SDW): 

J± = 2Ji (cos 5± + cos 5f) ± 2 J 2 (cos + Sfj + cos (<$J - Sfj + cos 25± + cos 28f) . 

(22) 

The quantities 8 yjZ are differences of the order parameter phase on neighboring chains. When 
J 2 is zero all 5's are zero, too. Otherwise, they have to be determined to maximize the 
effective coupling constants J±. Simple calculation shows that 5^ z deviate from zero when 
I J 2 | grows bigger then J1/6: 

r4J! + 8J 2 if J 2 > -Jl/2, r4J!-8J 2 ifJ 2 <Ji/2, , . 

+ I -4J 2 - J?/{3J 2 ) if J 2 < - Ji/2, " I 4 J 2 + J?/(3J 2 ) if J 2 > Ji/2. 1 ^ j 
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Of course, J\p > in our system. We will consider the case of negative J2 for completeness 
of the presentation. 

We did not put the products of the Klein factors in flT9l). This is permissible as long as 
the superconductivity and SDW do not co-exist. Then it is possible to use the following 
representation of the Klein factors: 

Wilt = i, VhVri = i, VliVri = i, Vr^Vli = i- (24) 

The averages of cosine of $ c and G c have to be found self-consistently. From the dimensional 
analysis: 

,a N l/(2/C c ) 

(cosV2^e c )«^J . (25) 
The required self-consistency condition reads: 



-a x 2-l/(2iC c ) 

J+ ( ~r~~ ) ocA s . (26) 



■s 

a; 

These two equations give the following expression for the gap: 

. (27) 

This is the equation for the gap due to the formation of the superconducting ground state. 
The equation for the SDW gap is: 

1 

A c oc A s (j=- J . (28) 

The zero-temperature transition between the superconductivity and SDW occurs when these 
two gaps are the same: 

1 1 

/ J,\ 2-l/Kc (J \2-Kc 

a '(ai) KA -fe) ' < 29) 

The critical temperature T c for a particular set of parameters equals to the biggest gap: 

T^A^U^,^) (30) 

The equation ( PU[ ) also gives the position of the tricritical point where the superconductivity, 
SDW and the phase without the charge gap co-exist. 

Now it is possible to construct the phase diagram. We will plot it on (J 2 , T) plane. 
It will be explained later how such diagram can be related to the (p, T) digram measured 
experimentally. One has to fix the value of /C c in the interval ( |15|) and r = J\/A s in the 
interval < r < 1. The value of J2 is constrained to: 

- Ji < J 2 < J x . (31) 
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The phase boundaries are determined by (|30D. The tricritical point is given by (p9|) . It 
is always located at right-hand side (J 2 > 0) part of the diagram. This is a consequence 
of (|15|). The point of T = transition is located right beneath the tricritical point. The 
diagram itself is presented on fig.l. The phase with no broken symmetry is denoted 'Spin 
gap'. 

Usually, experimental results are presented on the (p, T) plane. In order to establish 
connection between fig.l and the experiment we make three assumptions: (i) the coupling 
constants Ji i2 are decreasing functions of the transverse lattice constant bj_; (ii) J 2 is more 
sensitive to change in bj_ than J\ in the sense that J2/J1 is a decreasing function of bj_; 
(iii) both pressure and temperature affect J 12 through variation of b±; all other effects of 
temperature and pressure are unimportant. The supposition (i) immediately implies that 
J12 are increasing functions of the pressure and decreasing functions of the temperature. If 
we neglect the dependence of J\ on b± = b± (p, T) at all then (ii) is trivially satisfied. In 
such a case by applying pressure at zero temperature we move the system from SDW at low 
values of J 2 (ambient pressure) to TSC at higher values of J 2 (elevated pressure). Thus, the 
structure of the phase diagram at T = is the same on both (p,T) and (J2,T) planes. 

At non-zero temperature we have to remember that the couplings Ji j2 are not inde- 
pendent of temperature. Due to the thermal expansion their (bare) values decrease as the 
temperature grows. To compensate this effect of the thermal expansion extra pressure is 
required. Therefore, the tricritical point moves to the right. This effect is commonly referred 
to as re-entrance. It is a well-established feature of the experimental phase diagram of Q1D 
organic compounds ||. 

Fig.l suggests that the critical temperature of TSC grows with the pressure. This contra- 
dicts the experimental data for the organic superconductors Jl|] which show that T c decreases 
for p > Ptcp where prep is the pressure at the tricritical point. We will address this issue 
in Section VII. 

Let us study the validity of the mean field approach. We will evaluate the energy correc- 
tion introduced by the operator AH,. This correction has to be smaller then the mean field 
energy E^ oc LN±A 2 /v c (L is the length of chains, N± is the number of chains). The first 
order correction is zero: (AH) MF = 0. Second order contribution is given by the formula: 

£ (2) = L J dxdr (AH(0, 0)AH(x, r)) MF , (32) 
AH(x, t) = e-J^MF^l^^eT/WMF^'). (33) 

The problem at hand is reduced to the estimation of the bosonic correlation functions. When 
x/v c or r exceed A" 1 the integrand of (^) vanishes exponentially fast. Thus, the integration 
can be restricted to A^ 1 < \Jt 2 + x 2 /v 2 < A^ 1 . At such small distances the effect of the 
charge gap is negligible. This simplifies our task even further: the bosons may be considered 
to be free. The evaluation of integral ( |3"2"D gives: 

oc (L/v c ) ( £ ) £ J;r«- " I{ : oc LN ± ^, 
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(34) 
(35) 



_J0 forTSC 
q ~ Ul 1 (7T + 5j / ,7r + 4) for SDW { ' 

Two-dimensional vector Rj shows the position of ith chain. The quantity z can be considered 
as an effective co-ordination number of a chain. For models with the nearest neighbor 
interaction only and q = the value of z coincides with the co-ordination number of a 
chain. In the above expression for a constant independent of A c has been omitted. We 
see now that the condition for the mean field ground state to be a good approximation is 

z » 1. (37) 

This condition is satisfied if the number of transverse directions is big or if the transverse 
interaction is sufficiently long ranged. Inequality Q57| ) is consistent with findings of [TT 
where corrections for the inter-chain mean field theory were also evaluated. 



V. EFFECT OF LONG RANGE COULOMB INTERACTION 

The phase diagram on fig.l is a mean- field result. It ignores the contribution of the 
Goldstone mode to the free energy. Such contribution is unimportant at T = 0. At T > 
those modes get excited and increase the entropy of the system. Due to the entropy of these 
modes the first order phase transition line separating TSC and SDW phases bents to the 
right creating re-entrance region even on (J 2 / ' J\,T) diagram. 

The reason behind the re-entrant behavior is difference between the Goldstone mode 
contribution to the entropy of TSC and SDW phases. The long-range Coulomb interaction 
opens a gap in the spectrum of the Bogoliubov mode (Goldstone mode in the supercon- 
ducting phase) but leaves some SDW sliding modes without such gap. Therefore, even 
if at a given value of Ji 2 the ground state energy of the superconducting phase is bigger 
then that of SDW the entropic contribution of the sliding mode may trigger the transition 
into SDW phase as the temperature grows. On (p, T) phase diagram both thermal expan- 
sion discussed in the previous section and the Goldstone modes work together forming the 
re-entrance region. 

To study the effect of the long-range Coulomb interaction the following term has to be 
added to the total hamiltonian: 

Hr - M£ V fdxdx' ^'M^M fa v Or f ^ ^ ^ j, «t (38) 

Here gc oc e 2 /(v c e) is the interaction constant. Simple power counting suggests that (|38|) is 
a marginal operator. Provided that gc <C 1 the effect of the long-range Coulomb repulsion 



on the mean-field gap (27) and (28) is negligible. However, it changes the dynamics of the 
Goldstone modes significantly. In SDW phase the long-wavelength hamiltonian of the sliding 
mode is: 

trsdw _ f dk\\ (Pk± 2 dw 

H -J^j2^ bA ' (39) 
#k dw = ~ U/C c H k Ht + L/c-i*} + ^j-ki + ^gcbl'j^- | ^ | j (40 ) 
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n k = fc|,e k , (4i) 



s 2 dw cx J_A s 6i(cosV2^$> cxA^i. (42) 



Here we used an approximation cos $ ~ — ($ 2 /2) (cos $) to simplify the transversal coupling. 
The meaning of the constant f S( j w is the velocity of the sliding mode in the direction normal 
to the chains. From this expression the dispersion of the Goldstone bosons is found to be: 

k 2 

u 2 = v 2 k 2 + v 2 dw k 2 ± + 2IC c v 2 g c bl 2 1 ^L_ (43) 
The structure of this formula becomes more clear if we re-parameterize it: 



^2 = / vj + v 2 sAw + vj v 2 dw 2(j) \ k 2 + 2 JC c v 2 c g c bZ 2 cos 2 0, (11 ) 



k 2 = k\ + k\ , tan0 = k±/k\\. (45) 

We see that the gap in the sliding mode spectrum closes at = tt/2, that is, for mode 
momenta normal to the chain direction. 

Similar procedure for the Bogoliubov mode gives: 

u> 2 = + cos 2^ (k 2 + 2IC c g c b- ± 2 ) , (46) 

v 2 c oc J+A B b 2 ± (cos v^e) 2 oc A 2 c b 2 ± . (47) 

Unlike SDW sliding mode, the Bogoliubov mode has a gap for any direction of its momentum. 
The difference stems from the fact that for TSC case the long-range Coulomb interaction is 
expressed not in terms of the field variable but rather in terms of the conjugate momentum. 

1 /2 

The gap is lowest for <p = tt/2. Its value is A c g c . 



VI. SO(4) TRANSFORMATION 

In this section we will show that at some value of the model parameters the symmetry 
group of the hamiltonian (|IB|) is enhanced. One can observe straightforwardly that ( |T3"D is 
invariant under the global shifts 

Qci(x) -> e ci (x) + c, (48) 
$ ci (x) -> $ ci (x) + c. 

The first of these equations corresponds to the gauge transformation, the second corresponds 
to the translation along the chains. Thus, the symmetry group of the hamiltonian is U(l) ® 
U(l). We now show that if the additional conditions: 

K c = 1, jff W = i 1 )' ! .I;y (49) 

are imposed (where i and j are sites in the 2D bipartite lattice formed by the chains and 
^—\y-j — i jf i anc i j are on the same sublattice and —1 if on opposite sublattices) then the 
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hamiltonian symmetry group is SO (4). To demonstrate this we construct an explicit set of 
SO (4) transformations which leave the hamiltonian invariant. As a first step it is convenient 
to re-fermionize the hamiltonian using (|5|) with a c instead of a. We do this first for the 
inter-chain interaction term: 

n eS± = SvrV ]T (VI^R^tAb- + Ai^HiAu^u) ( 5 °) 

-Ji- (4 t A^^Ti + ^uAiAu)} e-^^ + h.c. 
To show the SO (4) structure of this operator we define a four component spinor 

*t = (<i/t = (^ lt ^ ^ (51) 

and 2x2 matrices: 

9 = ^1=1^ ^40- (52) 

With this notation the hamiltonian density 7i e tt± can be re-written as: 
H eS± = 47r 2 a c £ ((^ w + ^) tr&r^K + (j* dw - J*/) tr^j) e- 1 ^*-"^) + h.c. (53) 

We now define a set of transformations which (i) leaves the spin boson part of (|53|) invariant 
and (ii) acts as SO (4) on matrices g. To do this we introduce matrices which act on the 
four-component spinor (|51"D. These matrices are determined by: 

<V = (-i/2)[7M,7»], (54) 

' CTfc \ / (Ji 



^ = E e «*( ff J. *« = U -<J" (55) 
with 7- matrices of Dirac theory Jl2| : 

7<=(£ ,1 = 1,2,3, To =(J J), 75 = 7o7i7273 = ("J \) ■ (56) 
Consider the following transformation group of the spinor \l/: 

= exp (i ^ ^V<v) * (57) 



with tu^jy being c-numbers. The form of ( p4| ) and flog ) shows that \&l and ^r transform 
independently under (|57|): 

*L = (58) 

*k = (59) 



exp t<V<V J = y- ) > ^> ^ ~ unitary. 



(60) 
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Because the two L components of ^ have the same spin boson part (i.e. same factors of 



exp(±i-y/7r/2 ( l ) s ) and exp(±iy 7r/20 s )) as do the two R components and the transformation 
( j57|) does not mix L and R, the spin boson is left invariant under (|57|). Therefore, spin boson 
operators commute with this transformation, as (i) claims. 

It is routinely proven in introductory courses on the quantum electrodynamics that the 
object = \E r '7^ 1 I r transforms as a four- dimensional vector under (pTD . Therefore, the 
transformations ([57]) comprise SO (4) group. Unlike the theory of Dirac equation we do not 
need to deal with R^. To prove (ii) it is required to know the action of these transformation 
on g. As a consequence of (^) and fl59|) this action is given by: 

9 = UgVK (61) 

On a bipartite lattice a global SO (4) rotation can be defined by specifying matrices U\ and 
Va on sublattice A and assigning Ub and Vb according to: 

U B = cr z U A cr z , V B = (t z V A (t z . (62) 



If equation (f49|) is obeyed then in ( |53|) the term proportional to ( J^ w + J| c ) is not vanishing 
only for i and j on different sublattices and the term proportional to (Jfj — Jf?) only for i 
and j on the same sublattice. So the total expression fl53|) is invariant under (0). 
The same re-fermionization procedure applied to 7i c gives: 

n c = £ ^ (£c + l/Kc) v c J2 (CAVw - ^LAVw) (63) 

i / cr 

-4tt (/C c - 1//C C ) (n LTi + n Lil ) (n RT , + n Rll ) - {K c + 1//C C ) Vs ((V$ s ,) 2 + (VO si ; 



£ i-^ (/C c + 1//C C ) Wc (^75^^) + 4tt (/Cc - l/JQ tig t a z gja z 



2 

(/C c + 1//C c ) v s ((V$ si )* + (VG 



2 , ^T7A 

si ) 



The first and the last terms are invariant under SO (4) rotations. The second term is not 



invariant, but vanishes if (|49 ) holds. 

To summarize, for generic parameter values the symmetry group of the hamiltonian 
(|T3|) is U(1)®U(1) but if extra conditions ( f49D apply then the symmetry of the low energy 
hamiltonian is SO (4). The symmetry transformations connect smoothly the superconducting 
ground state and the SDW ground state. A similar situation takes place in U < Hubbard 



model on a square lattice 0] and in a certain Q1D system near charge density wave-singlet 
superconductivity transition |TT|. In the latter paper the symmetry was investigated using 
non-Abelian bosonization technique and found to be SU(2). On phenomenological level 
our model is an analog of SO(5) theory of high-T c superconductors proposed by S.-C. Zhang 



14]. The long-range Coulomb interaction destroys this symmetry (as seen, for example, 



from (|38|) ). 



VII. DISCUSSION 



In this paper we proposed a many-body theory of the transition from the spin-density 
wave to the superconducting state in Q1D materials. The essential component of the pro- 
posed mechanism is the existence of the high energy spin gap on a single chain. This gap 
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is a many-body phenomena unaccessible from the mean field theory. At the temperature 
lower then A s all intra-chain interactions and particle hopping are modified by the presence 
of the spin gap. Only the exchange interaction and the Josephson tunneling survive. The 
competition of these two determines the low-temperature phase of the system. Our picture 
allows us to produce a phase diagram which is qualitatively similar to that of the organic 
superconductors. In particular, we were able to explain the re-entrant region of the phase 
diagram. 

In our approach the superconducting phase is deeply connected to SDW. This connection 
manifests itself through the existence of the quantum symmetry between the superconduct- 
ing ground state and SDW ground state. In any real system this symmetry is violated. 
Depending on the violation we have either TSC or SDW as a ground state. As we demon- 
strated, the physical long-range Coulomb interaction further damages the symmetry by 
modifying the Goldstone spectrum of TSC and SDW in different ways. If this symmetry 
can be observed experimentally is unclear. 

It is also important to note that the superconducting phase is stabilized in a system with 
purely repulsive interactions. The total energy is reduced because the energy of transverse 
hopping is smaller in the superconducting phase. This is in obvious contrast with BCS 
theory where the optimization of total energy is achieved by lowering the potential energy 
of the electron-electron interaction. Even more interesting, in BCS superconductivity the 
pairing operator is only marginally relevant while in our system it is always relevant. Thus, 
T C BCS oc exp(— 1/g) where g is the coupling constant. In our case g oc J+/A s and T c is 
proportional to some power of g. In this sense, the proposed mechanism is 'high-temperature' 
superconductivity. 

A crucial prerequisite for experimental implementation of the mechanism is a strong 
Ising anisotropy of the spin-spin interaction. Although, in the organic superconductors this 
interaction does posses some anisotropy, the latter seems to be quite weak |15|] . The situation 
becomes even more aggravated if one realize that the size of the spin gap (|D|) is exponentially 
small at weak anisotropy. There are two possible ways out of this problem. First, one can 
try to improve the estimate Second, one may consider the system with the small value 
of A s . We believe, that the latter is a correct route toward the realistic description of the 
Q1D organic metals. An indication that the regime with small A s is experimentally relevant 
comes from dependence T c = T c (p). Experiments show that the critical temperature is 
a decreasing function of the pressure. The situation can be explained qualitatively if one 
assume that at high pressure the system becomes 3D anisotropic Fermi liquid. Since Fermi 
liquid with repulsive interaction is stable against transition into the superconducting phase 
the superconducting properties deteriorate under pressure. The spin gap also has to close 
before the system can become Fermi liquid. Thus, the compound inevitably enters the 
regime where the spin gap is comparable or smaller then the transversal interactions and 
the transversal hopping. Capturing this regime is a significant theoretical challenge. 

The major reason why the system with small A s is difficult to describe is the presence 
of the unquenched operator of the single-electron transversal hopping. When expressed in 
terms of the bosonic fields SjC and $ SiC this operator has distinct non-local structure: it 
creates kinks in both O s c and $ SiC fields. This makes it impossible to apply directly our 
method to such systems. However, an important conclusion can be drawn. Q1D metal 
without well-developed spin gap on every chain is close to an array of transversally coupled 
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LL. This puts us in a situation studied before with 'g-ology' approach ||. We already 
mentioned in Introduction about its failure to describe the superconducting phase. Our 
study suggests that the breaking of the spin-rotational invariance can help recover this phase. 
Indeed, the scaling dimension of the spin-density and charge-density wave susceptibilities is 
ddw = /C c + /C s . That of the Cooper pair susceptibility (longitudinal triplet and singlet) is 
d sc = 1//C C + /C S . In case of full rotational invariance (/C s = 1) only one of these two is smaller 
then 2. However, if /C s < 1 both can be less then 2 simultaneously. In this situation there 
are several relevant operators in the system. As we have seen, at the ambient pressure the 
ground state is likely to be a density wave since the exchange interaction is more relevant 
then the Josephson coupling operator (e dw < e sc if /C c < 1). At the elevated pressure the 
superconducting state benefits from next-to-nearest-neighbor processes and may become the 
ground state. 

Despite the appeal of this discussion we have to issue a warning. The arguments we made 
in the previous passage are based on the knowledge of LL scaling dimensions of different 
symmetry-breaking operators. In the presence of the unquenched single-electron hopping 
those dimensions may change in some unknown manner. That is why in this paper we chose 
to study the case of well-developed spin gap where the single-electron hopping is irrelevant. 

To conclude, we propose a mechanism of SDW and superconductivity in Q1D mate- 
rials. This mechanism stabilizes the triplet superconducting state without any attraction 
between the electrons. The effective hamiltonian possesses approximate SO (4) symmetry 
which connects TSC and SDW ground states. Depending on how this symmetry is broken 
the low-temperature phase is either TSC or SDW. 



VIII. APPENDIX 

In this Appendix we will derive the effective low-energy hamiltonian (|i~3D. First we 
recognize the fact that in the presence of the well defined spin gap A s the transverse single- 
electron hopping is ineffective: when a electron is added to or subtracted from a chain it 
creates a soliton in the field <3> s whose energy is at least A s . Yet, correlated hopping of two 
electrons can survive. For example, an addition of a Cooper pair to a chain does not induce 
a kink. This can be easily seen from (|14|): the operator A^o does not contain exponentials 
of © s which create kinks. The same is true for an addition of a particle-hole pair given by 
operator (0). In fact, these two possibilities are the most relevant inter-chain processes. 
Their competition dominates the low-energy properties of our system. 

To obtain the low-energy description we have to eliminate all states whose energy is 
higher then A s . A particular version of the elimination procedure we will use here is very 
similar to the Schrieffer-Wolf transformation known in the theory of Kondo effect. The idea 
is to split the Hilbert spaces of a single chain into two subspaces. The charge subspace W C j 
for the chain j is spanned by the vectors of the form \ip C j) |0 S j) where |0 SJ ) is the ground 
state of the spin boson. The effective hamiltonian density on the subspace W c = % W C j 
can be defined by its action on \ip) G W c : 



/ 



J2H t - H ± V^ £# sj + H cj - E V^H ± V c |V) =E\1>), (64) 



V 
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Here V c is the orthogonal projector on W c and V^r = 1 — V c . 



Unlike usual hamiltonian (34) depends on the eigenvalue E. We will demonstrate that 
this dependence is weak and for our purposes ([54]) can be treated as an ordinary hamiltonian. 
The expression ([64] ) is quite general. It is possible to simplify it by noting that the partial 



matrix element (0 SJ | H± |0 SJ ) is zero. This allows to drop the projectors from formula (|64]): 
/ 



V 



J 




^2H ci -(H ± [J2H s , + H c ,-E] HA \\ij c ) = E\ij c ), (65) 



where (. . .) s denotes the partial matrix element of an operator with respect to the spin boson 
ground state. We will calculate this matrix element below. 
It is convenient to re-write the inverse operator as follows: 



-i 



{H ± \J2H sj + H cj -Ej H^=J Q du^H ± e- u ^^^JH ± j (66) 

J "0 

where we took into account that, to the lowest order in t±, (j2j H c j — E^j \ijj c ) »s and as 
well as that (H±(0)H±(u)) s vanishes exponentially for u > A~ . The calculation of the 
partial matrix element (H± (0) H± {u)) s can be reduced to the calculation of the Green's 
function for the spin bosons: 

(H ± (0) H ± (u)) B = J dxdx' ]T (Hjaj (x, 0) H ±ij (x', u)) s (67) 
dxdx'J2 t2 ±ijJ2Jl 



/ s 

The operator (p pa is defined by (^). The first term in the brackets corresponds to a tunneling 
event in which an electron and a hole hop from chain i to chain j. The second term 
corresponds to tunneling of a Cooper pair (two electrons) from j to i. For every pair of 
chains i and j there are 64 terms in the expression (|67|) . Yet, only a fraction of them has 
non-zero partial matrix element (0 S | . . . |0 S ). For this matrix element to be non-zero the term 
must be 'neutral' with respect to kink creation operator exp (iyfjr/2Q a ). If at u = a kink 
(anti-kink) at ith chain is created then for the term to be 'neutral' an anti-kink (kink) has 
to be created at u > 0. This condition immediately reduces the number of possible terms. 
The surviving terms have the form (T^j where T is one of the following: 



V>l T (z,0)V> RT (a:',u), ^i l (x,Q)^ Kl (x , ,u); (68) 
^(x,0)ip Ri (x',u) 7 ^(x^O^l^x^u); (69) 
^ L1 (x,0)^ Ll (x',u), ^(x,0)^ Rl (x',u); (70) 
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■0LT ( x > °)^lt ( x ' > M ) > V'rt fa °) V>Rt fa, w ) ! 
V4,4.(a?» 0)Vt| fa> M )' Wufa 0)^ Ri (x, u). 



(71) 
(72) 



On this list only two first lines are of interest to us since they are the most relevant operators. 
Let us calculate the matrix element 



(73) 



V>L T ;fa 0)V> LTi (:r, 0)V4 Ti (x', w)^ RTi fa, w) 
^LTifa 0)V> Rti fa, u)) (V^fa 0)^ RTi (x', u) 



_J_ r) i -J „i^/^(e c^ (x,0)+^ cl (x,0)) -i^y^(0 c i(x\u)-'P c^ (x' ,u)) 



xe J 



i v / 7r/2(e cj (a;,0)-$ cj (a:,0)) c ,-i v /7r/2(e cj (a;',«)+<J' cj (a;',M)) 



|£ S (X - X, «) 



0, 



e i- v / 7 r/2(e si (a;,0)+*si(a;,0)) e -i- v /7r/2(e si (a;',u)-$«(a;','u)) 



0. 



(74) 



The correlation function Q s decays exponentially for \x — x'\/v s and \u\ above A" 1 . Unfor- 
tunately, the operator (|73| ) is quite complicated due to its non-local structure. An obvious 
way to simplify it is to approximate Q s by a delta-function. This step is acceptable for the 
low-momentum boson whose wave length is bigger then v c /A s . Yet, the presence of high- 
momentum bosons render this procedure unsatisfactory. This obstacle can be overcome by 
eliminating high-momentum bosons from the effective theory. The procedure itself is very 
similar to the elimination of the spin bosons: in addition to averaging over the spin bo- 
son ground state it is required to average over the ground states of high-momentum charge 
bosons. After this procedure ( |T3[ ) is reduced to: 



J „i ^ ^ 



dxdx'du ((v4, ti fa 0)ip hy (x, 0)^ RTi fa , u)ip m (x', u) 
VUV ^ VU J dxe ^(*M-^(*m J dx ' du \g s{x - x > iU) f\g> { 
g>(x,u) = (o> 



s,c> 
2 



ivV5(e> (*,())+»> (*,())) -iy/7/2(e>(x',u)-*>(x',u)) 



(75) 



(76) 



A ,K c /2 



oc ( -f ) ( A^ 2 + x^/vl 



(l//C c -)C c )/4 



In this formula the superscript '>' is used to denote the quantities associated with the high- 
momentum charge bosons. Putting the expression ( [75] ) back into ( |B"6"D and performing the 
integration over u and x — x' one gets an inter-chain interaction term of the form: 



jsdw^i I . . - 1 . 



dxe iV2^(<b ci (x,0)-$ cj (x,0)) 



JS dw = t 2 



An 2 a 2 



dxdu \Q s (x, u)Q^(x, u) 



(77) 
(78) 



_1 = A s /f c , r = \Jv? + x 2 1 1 



The evaluation of this integral requires knowledge of Q s at the intermediate values of r: 
A -1 < r < A" 1 . This knowledge is unavailable. To circumvent this problem we neglect any 
anomalous dimension to Q s and approximate it by (Ar) -1 / 2 . This gives: 
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(79) 



Other terms can be processed in a similar fashion. The resultant inter-chain interaction is: 
H^± = E {vUVkAvU + '/I.'/h.'/h.'//.) / dxa~ c l e^--^ (80) 

ij 

~J SC {vhvkAvU + liAWt) / dxa; 1 ^**-**) + h.c. 

The value of J sc can be obtained without any calculation if one takes the duality transfor- 
mation of the charge hamiltonian: C <-> $ c , /C c «-> 1//C C . The quantity dual to J sdw is J sc . 
A consequence of ([f|) and the duality is the following estimate: 

+2 / A \ (/Cc+l//C c )/2-l 

J- oc ^ f ^ cx J sdw (81) 



A s V A 



Thanks to the duality, we can make a more accurate statement about these two coupling 
constants: J sdw = F(/C C )J SC , where the function F(JC C ) is unity when K. c is unity. Away 
from this point J sdw /J sc = 0(1). 

The cut-off of this effective theory of the charge bosons is now equal to a^ 1 . 

The products of Klein factors do not commute with each other. Thus, in general, they 
cannot be dropped. Yet, within an inter-chain mean field approximation the contribution of 
the Josephson tunneling term in SDW phase is zero and, conversely, the contribution of the 
exchange term in the superconducting phase is also zero. In such a case, it is permissible to 
ignore Klein factors. 
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FIGURES 




FIG. 1. Mean- field phase diagram of our model. Solid lines mark the location of the sec- 
ond-order phase transition into SDW or TSC. Dash line corresponds to the first order phase 
transition between SDW and TSC. Label 'Spin gap' denotes the phase in which the charge bosons 
are disordered. The gap in the charge sector (equations ( |27| ) and (28)) are shown by dotted lines. 
The unphysical region J 2 < is shaded. 
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